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Abstract 

The problem of the vortex structure in the lattice model of superconductors has been reduced 
to the nonlinear map problem characteristic for the fractal theory. 
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I. INTRODUCTION 



The lattice magnetic translation symmetry of superconductors near the upper critical 
magnetic field was shown to be crucially important for the critical thermodynamic properties 
U, |2j, ||. Attepts to construct the incommensurate vortex structure were undertaken in 
0]. The small parameter (3 = $/$o (where $ = Ha 2 and $o — ff are respectively the 
magnetic flux through the lattice plaquette a 2 and the London quantum of magnetic flux) 
was essentially used there that predetermined very weak dependence of the sample free 
energy on the vortex configuration. Here I formulate this problem in a different way, so that 
smallness of (3 will be not assumed. 



II. BASIC EQUATION 

Confining myself here by the mean field approximation I consider the two-dimensional 
model. Apart from neglecting of fluctuations this means that an entanglement of vortices 
will not be accounted for in this paper. The free energy functional for the lattice model 
reads: 
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where J is the tunneling integral, the lattice site coordinates can be written as x = ma, y 



na, (j) nm stands for the complex conjugated order parameter, A is the vector potential, 
integration of it is carried out along the straight line connecting the sites nm and n'm', 
-. Differentiating eq. ([D with respect to <p nm we obtain the lattice Ginzburg- 
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Landau (LGL) equation for a superconductor at a strong magnetic field: 
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Choosing the Landau gauge A = e y Hx and considering the simple square lattice case within 
the tight-binding approximation (with J standing for the nearest-neighbor tunneling inte- 
gral) we can write eq. @ in the form 

J {4> m +i,n + 4>m-i,n + 4>m,n+i exp [-i2nam] + (f) m ,n-i exp [i2iram]} + 
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Eq. ([3D can be transformed to the simpler form. Notice that eq. (|j) is a nonlinear equation; 
therefore, the Fourier transforming is not particularly useful since it will convert this equation 
into an integral-difference one. However, the specific form of nonlinearity characteristic for 
the Ginzburg-Landau equation allows us to make a substitution that is usually done in the 
theory of the linear Harper equation |J: 

4>mn = u m exp (mm) . (4) 

Then eq. (|3]) takes the form 

u m+ i + w m _i + 2 cos (2nm(3 - k) u m + ru m + g \u m \ 2 u m = 0. (5) 

Eq. ^ can be named the nonlinear Harper equation (NHE) or the Harper-Ginzburg-Landau 
(HGL) equation. The rest of the paper will be devoted to a discussion of possible conse- 
quences of this equation. 

III. ANALYSIS OF THE NONLINEAR HARPER EQUATION 

Eq. [| is very peculiar one. It contains two mechanisms, either of which can lead to 
forming of irregular self-similar structures that have fractal properties. If we shall proceed 
in the spirit of the Abricosov vortex structure theory, we have firstly to linearize the equation. 

u m+i + w m _i + 2 cos (2iimj3 — k) u m + ru m = 0, (6) 

This equation differs from the linear Harper equation || in the same way as the Ginzburg- 
Landau one for a superconductor at the magnetic field differs from the Schroedinger equation 
for a charged particle at the magnetic field. Therefore, the line H C 2 is determined by the 
condition r = q, where e\ is the Harper multiband spectrum lower edge. The dimensional 
crossover takes place near the lowest band bottom at (3 rational 0]. The flux structure 
is completely pinned by the crystal in this case. A subsequent decrease of temperature 
or magnetic field, i.e. an increase of |r|, leads to appearance of fractal structures. The 
Harper operator spectrum is known to form a Cantor set at irrational (3 and, therefore, 
the dimensional crossover is absent in this case. Scaling properties of this operator have 
been studied using the renormalization group technique || 0. Exact solutions with a use 
of the Bethe Ansatz and quantum groups can be found in p. Mult if r act al properties of the 
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Harper operator eigenfunctions and eigenvalues are discussed basing on numeric calculations 
and on the Bethe Ansatz in the work by Wiegmann with co-authors ||. Interesting results 
have been also obtained in Unfortunately, the Bethe Ansatz rigorous approach is not 



extended to the nonlinear case; it is a challenge now for the mathematical physics community. 
An approximate solution of the linearized Harper equation can be presented as a chaotic 
distribution of the Wannier-type functions, which can be considered as nucleation centres 
of the superconducting state. The order parameter amplitude can be found perturbatively 
in this case. On the other hand, Eq. (||) can be written in the form of the two-dimensional 
nonlinear map on the complex field: 
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IV. CONCLUSION 



The main inference from the work is that the fractal flux structure is controlled by 
two-dimensional maps on the complex field enumerated by the whole real numbers. This 
means that the two-dimensional flux lattice in the crystal lattice can be completely specified 
establishing the nucleation centre distribution along one say leal axis. In opposite to the 
continuum model case, this distribution is not periodic in general. The following distribution 
can be taken first approximation 

(f)(x,y) = Cexp ( x (8) 
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where u n are determined by some proper discrete map taking fractal effects into account. 

Notice that this sum is reduced at the condition u m = into the elliptic theta function so 
that this formula can be considered as its natural generalization on special number fields 
similarly to introduction of the basis functions, for instance. Along with the discussed above 
map, one can try the Chirikov standard map, which is equivalent to the Frenkel-Kontorova 
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model that is widely used for description of the incommensurability effects: 

I n+1 = J n - 2gb sin ( 2 7Cy n /b) 



(9) 

Dn+l = Vn + Ct + In+l) 



or, if the interaction between vortices must be taken into account, the generalized standard 
map 

f / n+1 = / n -2^sin(27ry„/&) 



Vn+i = Vn + a ~ p ln(l - (3I n+1 ) 
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